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NOTE ON THE SIMPLE DIFFERENCE EQUATION. 

Bt J. H. M. Wbddbkbubn. 

Boole on his Calculus of Finite Differences derives the Euler-Maclaurin 
sum formxila from the expansion of the operator A~^ = (e* — 1)~^ in a 
series of positive powers of D so obtaining the well known result 

(1) A-^g(x) = fg{x)dx - Mx) + £ (- D" (^^+2) ! ^^ + «(^)' 

where «(x) is an arbitrary periodic function of period unity. This series 
is usually divergent and is besides only applicable to functions which 
possess derivatives of every order. 
If instead of expanding (e* — 1)~^ in powers of D we use the series 

^'^) gi* _ 1 B 2^ r L-D + 2nid D - 2mn\' 

we obtain another form of the solution of the equation A/(x) = g{x), 
namely 

g{t)dt + 2 Z J cos 2m{x - t)g{t)dt + w(x), 

the lower limits being arbitrary constants whose imaginary parts are all 
equal and the path of integration a line parallel to the real axis. 

The series (3) has this advantage over (1) that it is usually convergent, 
the conditions of convergence being similar to those of a Fourier series. 
It is however in general discontinuous; for, when the lower limits are all 
made equal to zero, then if 

a; = 0, /(O) = - k(0), 

^ a; s 1, f{x) = 0, 

X = 1, /(I) = MO), 

1 S x S 2, fix) = g{x - 1), 



n ^ X S n + 1, fix) = gix - 1) + gix — 2) + • • • gix — n). 

As the derivation of the series (3) is purely formal, it is necessary to 
verify directly that it does satisfy the equation A/(x) = gix). Forming 

82 



NOTE ON THE SIMPLE DIFFERENCE EQUATION. 83 

A/(x) we have 

ri » /.i+i 

g{t)di + 2 Z J^ cos 2Tm{x - t)g{t)dt 

= - h\g(x + 1) - g{x)] + ^[g{x + 1) + g{z)] = g{x), 
since the series 

g{t)dt + 2 Z I cos 2?m(x - t)g{t)dt 

is the value of the Fourier series for g{x) at the end of the range of integra- 
tion. Hence (3) gives a value of A~V(x) when g{x) is a continuous func- 
tion capable of being represented by a Fourier series. 

This solution is of course not essentially different from the sum 
formula of Poisson: it is also closely related to the Guichard integral 
from which it may be derived. 

2. As an illustration of the use of our form of the solution we shall 
now derive from it the well known asymptotic series for log r(x) which is 
a solution of the equation 
(4) A/(x) = log X. 

Forming the series for/(x), we obtain as a solution of (4) 

log id« + 2 Z J cos 2im(x - t) log tdt 

/ iM , 1 , ov^ rsin2im(x - <), .T 
= (x -^)logx -x-h 1 -f-2i:- [ 2^ logfj^ 



= (x-§)logx-x + H-2f:J' 



sin 2im(x — t) 

X 

sin 2Tm{x — t) 



00 /»* 



2imt 



-dt. 



If we now change the lower limit of integration from 1 to «, the value 
of the series is only changed by a periodic function of period unity, so that 



Uix) = (x - ^) logx - X -I- A; - 2 E f 

1 *Jx 



' sin 2m(}C — ,^ 
at 



= (x-i)logx-x-h^ + 2i:(2^,x- + 2i:J^ (2^,).,. dt, 
where A; is a constant, is a solution provided the series used are convergent. 



(2im)2f*~a;'r(2im)2' 
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When X is real and positive, 

I Ql _ Iv r" cos27m(x - ^f, f dt 1 ^ 

' ^ ' - 14- J. (27m)^<^ '^^ = 4- Jx (27m)^f^ " a; Vi 

while, ii X = i + rfi, t = a + Tji, (ri =^ 0), 

I '^ ' = r jf {2mnc^ + v')^ \v\r (27m)^' 

The series therefore converges and 

2 °° 1 

fiix) ~ A; + (a; - i) log X - X +^l) (a^)2 

so long as x does not lie on the negative half of the real axis. It is easily- 
shown by the \isual methods that k must equal | log 27r in order that 
fi{x) may coincide with log r(x). 

3. Another solution of the simple difference equation may be obtained 
by means of the integral equation 

git)dt = J^ mdt, 

which on differentiation leads to 

g{x) =f{z + l) -fix). 

As is suggested by the solution of the Fredholm equation, we consider 
first the homogeneous equation 

Xi+i /•<» 

u{t)dt = \J Kix, t)u(t)dt, 

where the kernel K(x, t) is defined by 

K{x, = 0, aSt <x 

= 1, x St ^x + 1 

= 0, x + 1 <t, 

a being a fixed constant. This kernel has a continuous spectnmi: for if 
in (8) we substitute 

(9) u(x) = e*^ (& 4= 2njrt) 

we get, after a simple reduction, 

(10) X = ^ 



e" - 1 
so that X may have any finite value. 
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The theory of integral equations therefore suggests the expansion of 
I g(t)dt in a series of the form 

(11) rgit)dt = Sa.e**". 

If an expansion of this form exists, a solution of (7) is 

(12) /(a:) =Sa,--q4r-i«^ 

provided this series is convergent. 

If g(x) can be expanded in a uniformly convergent series 

(13) g{x) = S&iC*^, 
then 6, = «,&,•, and we may set for f{x) 

(14) fix) =:2^s^e'^ + o,{x). 

When g(x) is a real function which satisfies Dirichlet's conditions, an 
expansion of the form (13) always exists. In particular we have* 

where the summation extends over all the roots of 

(16) <T + h tan a<T = 0, 

and — a < X < a. 

None of the roots of (16) make e'* — 1 = 0, and as <r approaches 
infinity, a approaches the value (2n + l)ir/2a, so that if o is an integer, 
e'* — 1 remains greater than some finite quantity. Hence the form of 
f{x) given by (14) is convergent and a solution of (7) is therefore given by 

^(^) = 2(1+^ 1/^^)^^+^ (sin 2a.-2a<r)(e-*-l) }j '^^(^)^^+"(^) 
hx r ,^.j^,„ <T r ,., sin(r(x-<-^) ,. , ,, 

= W+^ j_^g(Od^+S ^^ _ 3-^ 2a.r 1/(0 sink '^^+"(^)' 

where the second summation extends only over positive values of a: 

• See for instance Picard, Traits d' Analyse, 2° ed., t. 2, p. 192. 



